Direct Imaging of Spatially Modulated Superfluid Phases in Atomic Fermion Systems 
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It is proposed that the spatially modulated superfluid phase, or the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state could be observed in resonant Fermion atomic condensates which are 
realized recently. We examine optimal experimental setups to achieve it by solving Bogoliubov-de 
Gennes equation both for idealized one-dimensional and realistic three-dimensional cases. The spon- 
taneous modulation of this superfluid is shown to be directly imaged as the density profiles either 
by optical absorption or by Stern-Gerlach experiments. 
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A keen interest focuses on resonance Fermion conden- 
sations by using ultracold atomic 6 Li or 40 K clouds. It 
is now agreed that the BCS superfluid state has been 
achieved by various experimental groups Q, 0, 0, 0, 0, IE 
.7], where magnetic field Feshbach resonance is utilized 
to finely control the atomic interaction between strong 
repulsive and attractive regions. There are not only ac- 
tive experimental and theoretical 0, investigations on 
this BCS-BEC cross-over phenomenon, but also pioneer- 
ing works on exploring the BCS pairing state, including 
spectroscopic experiments to extract the energy 

gap value as functions of temperatures and the coupling 
constant 0,0. 

In the current experiments, the prepared Fermion 
numbers of two species, where the hyperfine states \F = 
%,ra F — — |) and \F — |,m F = — |) are used in the 

K case are equal to maximize the transition tem- 
perature. Here by intentionally making these numbers 
unequal, we could achieve a spatially modulated BCS 
state, namely so-called Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state 0E3- 

This interesting state has been much sought out for 
long time in vain in a superconductor. One of the main 
difficulty to stabilize it lies in the fact that applying 
field to create population difference of up and down elec- 
trons inevitably induces diamagnetic current which acts 
a Cooper pair breaker. An exceptional case is ferro- 
magnetic superconductors whose internal ferromagnetic 
molecular field does not induce it. Even under this favor- 
able situation any firm conclusion has not been reached 
for ternary compounds ErRli4B4 and HoMofiSgflfil Il7j. 
One of the difficulties to observe it is the absence of an 
experimental method to see the modulation. The FFLO 
state is also discussed in a new heavy compound CeCoIns 
|l8| or even in hig h energy community in connection with 
neutron stars 

The present neutral Fermion systems are ideal to pur- 
sue this realization of the FFLO state because: (1) There 
is no pair breaking mechanisms such as the diamagnetic 
current or impurities which weaken the stability of this 
state. (2) As is shown below, fine adjustment on the 
population difference of two species, say, up and down 



Fermions, and attractive interactions is possible for the 
system to be in most favorable and easily observable con- 
dition for the FFLO state. 

The purposes of the present paper are to analyze the 
most suitable condition to observe the FFLO state by 
maximally utilizing the limited experimental resources 
for the present neutral atomic vapors and also to provide 
the experimental signatures of the FFLO state within 
the currently available experimental techniques. The ar- 
rangement of this paper is following: After introduc- 
ing the Bogoliubov-de Gennes formalism which allows 
us to describe the spatially modulated phase, we pro- 
vide briefly known analytic solution for the FFLO state 
of the quasi-one-dimcnsional system which is useful 
to grasp the general perspective of the present problem. 
Then, we go on showing a numerical calculation for the 
more realistic situation where the atomic cloud is con- 
fined by a three-dimensional harmonic trap. 

We start out with the pairing Hamiltonian 



H= dr 
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with g = 



4mH a 



being the attractive interaction. a(< 0) is 



the s-wave scattering length. An axis-symmetric trap po- 
tential V(r) = im(^r 2 +cu^z 2 ) is assumed. Here, to pre- 
pare the unequal population in two species a =t, I, the 
following chemical potential is used: \x a = fi + (o'3) ITi(T 5/x 
with the Pauli matrix 03 . The Fermion operators of the 
two species are denoted by V>£(r) an d VvM- The Bo- 
goliubov transformation with the mean-field approxima- 
tion is employed to diagonalize the Hamiltonian, that is, 

= EqK'Jq,! - "q^qj aild V 5 ! = EqK»7q,T ~ "q^qj 

within the creation and annihilation operators rjq ^ and 
77! CT . To this end, one can lead to the Bogoliubov-de 
Gennes (BdG) equation: 
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where K a = -|^V 2 + V(r) + gp a (r) with the self- 
consistent equation 

A(r)= flefi ^ Wq (r)<(r)/( £q ). (3) 
q 

The order parameter A(r) = ^ e fF (^i ( r ) ( r ) ) an< ^ 
the particle density of each component p-f(r) = 
E q l«q(r)| 2 /(£q), and p;(r) = £ q Mr)| 2 [l - /( £q )] 
where /(e q ) = \j(e e ^l kBT + 1) is the Fermi-distribution 
function. Since the finite shift of the chemical poten- 
tial (T3 5 p causes a distinction of the mean- field poten- 
tial for two species, the particle-hole symmetry is broken 
|2*o) . Therefore, the sum in Eq. J3J), is done for all the 
eigenstates with both positive and negative eigenener- 
gies. When calculating A(r), we have used a regularized 
coupling constant g c g to avoid the ultraviolet divergence 

In the one spatial dimension without a trap, which is 
approximate for a long cigar-shaped gas bJ z juj r — > 0, one 
admits an exact analytical solution Ijj, which has the 
form 

A(z) = Aifcisn(Aiz, kx). (4) 

Here, sn(z, k) is the Jacobi elliptic function with the mod- 
ulus k {k\ complementary modulus). Ai is the order 
parameter to be determined self-consistently. According 
to this solution at zero-temperature, the FFLO state is 
energetically stable for the relative population difference 



n 7r ep ' 

where the particle number of each component n a — 
J dvp a {v) and the total particle number n = J2 a n a . Ef 
is the Fermi energy and An is the amplitude of the order 
parameter of a uniform BCS state at zero-temperature. 
Beyond the above critical population difference the uni- 
form BCS state changes into the modulated FFLO state. 
Since Aq/ep = 0.2 ~ 0.4 in the present experiments, 
the critical population difference should be of an order of 
10% ~ 20% difference. 

This modulated phase accompanies the spin variation 
with half of the fundamental modulation periodicity. The 
excess density of the up-spin Fermions (n-f > njj periodi- 
cally accumulates at the zeros of A(z). The accumulation 
becomes sharp in space when either approaching the crit- 
ical population difference or stronger coupling region. It 
is illustrated in Figs. 11 and 12 of Ref. |lj| where A{z) 
and the associated spin modulation are seen for two typ- 
ical Sn values, one is near the critical boundary and the 
other far from it. 

Since each atomic species can be probed optically be- 
cause of different absorption frequencies, we can image 
each density separately in space. Thus it is crucial to ap- 
propriately choose the two parameters; \Sn\ = njj/n 



and Ao/ef in order to obtain a sharp absorption image. 
Another way to detect this spin modulation is to per- 
form Stern-Gerlach experiments, revealing the different 
particle density distributions for spin-up and spin-down 
Fermions. 

As for the modulation periodicity A relative to the 
mean particle distance L — 1/kp, 

A _ eg kK(k) 

L An 7T 

where K(k) is the complete elliptic integral of the first 
kind. As a rough estimate we find A = 3/im ~ 30/im for 
A/sp = 0.1 ~ 0.01 and n ~ 10 20 ur 3 . 

After having finished the order of magnitude estimates 
for various quantities based on the analytic solution, we 
present a more realistic calculation by numerically solv- 
ing the BdG equation (J2J) self-consistently for help to de- 
sign an experiment: We consider the axis-symmetric har- 
monic trap potential with uj r /2ir — 1.67cj z /27t = 1000Hz 
and the 6 Li atoms with the chemical potential p = 
I2.5ft,uj r , corresponding to n ~ 1100 atoms. Through- 
out this paper, A /sp — 0.35 and ksT/?uv r = 0.05 are 
fixed. 

The results are displayed in Fig. ^ where the spatial 
variation of the pairing field A(r, z) and the total density 
profile p(r,z) are shown. A(r, z) changes its sign at the 
plane z = and also near the edges of the pairing field 
(See the red curve in Fig. dc)). Since the system has 
an elongated ellipsoidal shape, A(r, z) is never a simple 
sinusoidal form. This is quite different from the above 
one-dimensional case where the total density is assumed 
to be uniform. In the present situation the total density 
is non- uniform because of the confinement. The order 
parameter is inevitably inhomogeneous intrinsically. 

As is seen from Fig. ^b) the total density along the 
z-axis is similar to that expected for Thomas-Fermi ap- 
proximation for the normal Fermions except for a small 
feature, a dip at z = 0. This is related to the local sup- 
pression of A(r). 

In Fig. Die) we show the cross-section A(r = 0, z) at 
the r = axis. As 5n becomes large, the modulation 
period tends to become short, thus the extra zeros appear 
near the edges. It is also noted from this that for Sn = 
0.32 the modulation period becomes shorter towards the 
edges where the local density tends to be dilute. 

In Fig. |2 we display our most important results, which 
are directly probed experimentally. The local magnetiza- 
tion m(r) defined by m(r) = pi(r) — /Qj,(r), which is the 
local population difference of the two species, is shown in 
Fig. 2(a). It is seen from this that m(r) exhibits peaks 
at the central plane z — 0, where a prominent peak is 
observed at r = 0, and also a small peak at the edge. 
Since we are treating a three dimensional object with 
an almost spherical shape, it is essential for experiments 
that the physical quantities integrated along the radial 
direction, namely the columnar densities should exhibit 
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FIG. 1: (a) The spatial profiles of the pairing field A(r, z) 
and (b) the total density p(r, z) for the population difference 
Sn = 0.15 and Ao/£f = 0.35. (c) The cross-section of A(r, z) 
at r — is shown for various population differences Sn. 
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FIG. 2: (a) Stereographic view of the magnetization m(r,z), 
and (b) its columnar density profile along the z-axis for 
Sn = 0.15. (c) The columnar magnetization profiles m(z) 
for various Sn — 0.07, 0.15, 0.32. The arrows indicate the po- 
sitions of the zeros of A(r = 0, z), which is displayed in Fig.0 
The number of the zeros is 1 (Sn = 0.07), 3 (Sn = 0.15), and 
5 (Sn = 0.32). 



a prominent feature to probe it. In Fig. ffi h) we dis- 
play the columnar density p a (z) = 2tt J drrp a (r, z) and 
m(z) = pi(z) — Pi(z), both of which are directly observ- 
able. Because the two species have different magnetic 
moments we can spatially separate the two species by 
passing it through a field gradient. In this Stern-Gerlach 
experiment, which is performed in spinor BEC 22], p-j (z) 
and Pi(z) in Fig. 2(b) are directly imaged. Another way 
to reveal these modulations is to utilize the fact that 
each species has different clock transition frequency. By 
taking an absorption image with a particular frequency 
tuned we can selectively probe these modulations. 

It is seen from Fig. 2(c) that depending on the popu- 
lation difference under a fixed A /ep, the magnetization 
m(z) shows quite different spatial structures: As 6n in- 
creases, the single peak structure for Sn = 0.07 changes 
into multiple peak ones, corresponding to the increment 
of the zeros of A(r, z) (also see Fig.dc)). 

It is instructive to see the gap structure in the FFLO 
state because the Fermionic excitation spectrum is dis- 
tinctively altered from the uniform BCS state. In Fig. |3| 
the local density of states (LDOS) for up and down spins 



defined by 

AMr,£H]>> q (r)| 2 ^-£q), 
q 

7V i (r,S) = ^| Uq (r)| 2 ( 5( J B + £q ), (7) 
q 

are shown. It is seen from Fig. 3 that: (i) Away from the 
center at r = and z = 0, the energy gap becomes nar- 
row as a function of z, corresponding to the low-energy 
quasiparticle excitations bounded near the surface of the 
cloud [2i|]. (ii) The localized LDOS appears near the 
Fermi level at z = and also near the edges of the cloud 
z ~ 8pm. In the spin-up (spin-down) case of Fig. |3f a) 
((b)) the large mid-gap state at z — situating below 
(above) the Fermi level is filled (emptied), meaning physi- 
cally that the excess spin-up Fermions are accommodated 
by keeping the overall bulk gap structure intact, (iii) 
These mid-gap states at z = correspond to the posi- 
tions of the zeros of A(r, z) shown in Fig.l. 

We can observe these mid-gap states by using stimu- 
lated Raman transition as an extra satellite. This tech- 



4 



N t (r=0,z,E) [arb. unit] (a) 




N+(r=0,z,E)[arb. unit] (b) 




FIG. 3: (a) The local density of states profiles iVf (r — 0, z, E) 
and (b) N±(r — 0, z, E) at r = for the population difference 
Sn — 0.15 and Ao/ef = 0.35. The large empty region in the 
center in both figures corresponds to the energy gap. 

nique is employed by Chin et al. 0], wn0 identify the 
energy gap in Fermi condensates. Note that this spec- 
troscopy is most useful for a cigar shape cloud to separate 
out the mid-gap and surface excitations. 

Finally we should mention a work by Combescot [24| 
who examines the stability problem of neutral atomic 
Fcrmions with different Fermi radii under the absent of 
a trap. Yet another superfluid phase of an anisotropic 
gap is predicted. The precise boundary between FFLO 
and this phase remains to be worked out in a realistic 
situation. 

We also mention that the present spatially modulated 
in FFLO shares a common physics of the so-called stripes 
paradigm, which includes spin density waves in Cr |25| . 
spin-Peierls system 26], and high-T c cuprates [2t|]. In 
these systems the order parameter of the staggered mo- 
ment is commonly characterized by the spatial modula- 
tion with sign change. The node of the order parameter, 
namely the domain wall feature corresponds to stripes in 
the latter. 



In conclusion, we have proposed an experimental way 
on resonance Fermion superfluid systems with unequal 
mixtures of two species to achieve the FFLO state. One 
attractive point to accomplish it is that under a fixed 
population difference, say 10% by moving out from the 
Feshbach resonance point, we can cross the boundary 
from the uniform Cooper paired state into the FFLO 
state. This boundary region is most favorable to im- 
age the modulated pattern since the magnetization has 
a sharp peak at the zero of the pairing field. The di- 
rect imaging can be performed either by species-selective 
absorption, or by Stern-Gerlach separation. 

One of the authors (K.M.) thanks W. Ketterle for en- 
couragement. 
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